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Abstract 

We study the statistics of the relative separation between two fluid 
particles in a random flow. We confine ourselves to the Batchelor regime, 
i.e. we only examine the evolution of distances smaller than the smallest 
active scale of the flow, where the latter is spatially smooth. The La- 
grangian strain is assumed as given in its statistics and is modelled by a 
telegraph noise. This is a stationary random Markov process, which can 
only take two values with known transition probabilities. The presence 
of two independent parameters (intensity of velocity gradient and flow 
correlation time) allows the definition of their ratio as the Kubo num- 
ber, whose infinitesimal and infinite limits describe the delta-correlated 
and quasi-deterministic cases, respectively. However, the simplicity of the 
model enables us to write closed equations for the interparticle distance 
in the presence of a finite-correlated, i.e. coloured, noise. 
In ID, the flow is locally compressible in every single realization, but the 
average 'volume' must keep finite. This provides us with a mathematical 
constraint, which physically reflects in the fact that, in the Lagrangian 
frame, particles spend longer time in contracting regions than in expand- 
ing ones. Imposing this condition consistently, we are able to find analyti- 
cally the long-time growth rate of the interparticle-distance moments and, 
consequently, the senior Lyapunov exponent, which coherently turns out 
to be negative. Analysing the large-deviation form of the joint probability 
distribution, we also show the exact expression of the Cramer function, 
which happens to satisfy the well-known fluctuation relation despite the 
time irreversibility of the strain statistics. 

The 2D incompressible case is also studied. After showing a simple gen- 
eralization of the ID situation, we concentrate ourselves on the general 
isotropic case: the evolution of the linear and quadratic components is 
analysed thoroughly, while for higher moments, due to high computational 
cost, we focus on a restricted, though exact, dynamics. As a result, we 
obtain the moment asymptotic growth rates and the Lyapunov exponent 
(positive) in the two above-mentioned limits, together with the leading 
corrections. The quasi-deterministic limit turns out to be singular, while 
a perfect agreement is found with the already-known delta-correlated case. 



1 Introduction 



Understanding the behaviour of the relative separation between two fluid par- 
ticles in a generic flow is a difficult task of paramount importance. The most 
intuitive application of this concept is represented by the passive scalar problem, 
where the Eulerian description in terms of fields has its Lagrangian counterpart 
just in the study of fluid trajectories pQ. 

While in the one-dimensional situation general considerations and results can be 
carried out (at least for smooth flows), only few cases are known to be solvable 
in higher dimension [5J [5] ; namely, a short-correlated strain (fully solvable) , a 
2D slow strain 0] and the large-dimensionality case [5]. 

In this paper we describe another situation which allows exact analytical 
calculations: the telegraph-noise model for the velocity gradient in the Batchelor 
regime 6 . This means that the relative separation between two fluid particles, 
R(t), evolves according to 

d t R{t) = a(t)R(t) , (1) 

and the Lagrangian strain (scalar or matrix) a(t) is assumed as properly made 
up of a telegraph noise [7]. We shall study both the ID (compressible) and the 
incompressible 2D cases. 

The telegraph noise is a stationary random process, a(t), which satisfies the 
Markov property and only takes two values, a\ and a 2 . The probability, per 
unit time, of passing from the latter state to the former (or viceversa) is given 
by vi (y2i respectively). In what follows, we shall consider some special cases, 
for which simple formulas hold. 

(i) If a\ = — ci2, then the process itself is stochastic but its square is determin- 
istic, keeping the constant value a 2 (t) = a 2 (with a = |ai| = | a.2 1 ) - We shall 
denote such processes with a star, e.g. a*(t). 

(ii) If the average (a) = {v\a\ + V20.2)/v vanishes (with v = v\ + ^2), then the 
autocorrelation takes the form {a(t)a(t')) — ^ v \ t ^ t \ {i/ l a\ + V2a 2 ,)/v. For any 
analytical functional F[a] (function of time), a simple "formula of differentia- 
tion" [7] then holds in this case: 

d t (a(t)F[a}) = (a(t)d t F[a\) - u{a(t)F[a\) . (2) 

We shall denote such processes with a tilde, e.g. a{t). Note that it is always 
possible to reduce to this case from a generic telegraph noise (in particular from 
one a*(t)), simply by subtracting from the latter its mean value. 

(iii) If both the properties in (i) and (ii) hold simultaneously, the maximum 
level of simplification is reached. For the sake of simplicity, we shall denote such 
processes not with a*(t) but rather with £(i). 

This type of noise has often been used to mimic more general colored noise 
in stochastic theories and has also been applied, e.g., to optics (see [8] and 
references therein). An interesting application of the "telegraph concept" to 
turbulence can be found in [9]. 

The paper is organized as follows: in section [2] we will study the ID case. 
In section [3] we will move to 2D, first describing a simple peculiar situation 
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(the hyperbolic flow, 13. ip and then the general incompressible isotropic case; 
the latter will be analysed at first in its complete form, then (respectively in 
subsections 13.2.11 and 13 . 2 . 2")) by means of a simplified dynamics and eventually in 
the quasi-delta-correlated limit. Conclusions follow in section^) The appendices 
are devoted to the quasi-deterministic limit ((TfJ and to some calculation details 

2 ID case 

The simplest meaningful model, making use of the telegraph noise, for the ID 
case, consists in assuming a strain function with the property (i), i.e. 

a(t) = a*(t). (3) 

Indeed, the problem is intimately compressible locally, i.e. within every single 
single realization of the noise, but we want the average (ID) "volume" not to 
explode or shrink asymptotically in time. If we denote the long-time evolution 
of the moments of the interparticle distance witt0 

(R n (t)) t_ £°° e^">* , (4) 

this means that we have to impose the constraint 

7 (1) =0. (5) 

From ([1]) and ([3j) , it is straightforward to obtain the system of the two first-order 
ordinary differential equations which describe the time evolution of moments. 
It is indeed sufficient to multiply the starting equation by the noise itself (taken 
at the same time) in order to get a closed system, thanks to the simplicity of 
the telegraph form. Namely, if we define «o = ( a *(t)) and rewrite accordingly 

a*(t) = a(t) + a , (6) 

we obtain 

d t R n (t) = nR^i^dtRit) = na*(t)R n (t) = n[a(t) + a ]R n (t) (7) 
and thus 

d t (R n (t)} = n(a(t)R n (t)) + na (R n (t)) . (8) 

We must then look for the evolution equation of (a(t)R n (t)), which is readily 
done by exploiting @: 

d t {a(t)R n (t)) = (a(t)d t R n (t)) - v{&{t)R n {t)) . (9) 

x It is easy to show that the dynamics defined by {TJ is such that R(t) always keeps the 
same sign. Therefore, throughout this section, negative separations can also be considered, 
but they are to be dealt with either in absolute value or as the ratio R(t)/R(0). 
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Taking (J7J) into account, the first term on the right-hand side becomes 

(a(t)d t R n (t)) = n(& 2 (t)R n (t)) +na (&(t)R n (t)) 

= n{[a*(t) - a ] 2 R n (t)) + na {a{t)R n (t)) 

= na 2 (R n (t)) + na 2 (R n (t)) - 2a {[a{t) + a ]R n (t)) 

+na {a(t)R n (t)) 
= na 2 (R n (t)) - na 2 (R n (t)) - na {a(t)R n (t)) , 

where we made use of the property stated in (i). Consequently, © rewrites as 

d t (&(t)R n (t)} = n(a 2 - a 2 )(R n (t)) - (v + na ){a{t)R n (t)) . (10) 



Equations (J8J) and (|T0|l constitute the system we were looking for. It is easy to 
recast the latter as a second-order differential equation for the quantity (R n (t)) 
alone, which can then be solved exactly: 

d 2 (R n (t)) + vd t {R n (t)) - n(na 2 + va a ){R n {t)) = 

K (/£"(£)) = } )nC -Wu' A +'in A a A +4nva +v)t/2 _|_ ^gCv r v l : +4n 2 ' d i +<±nvu a -v)t /2 

=*> 7(71) = (y/v 2 + 4ri 2 a 2 + 4nua - u)/2 (11) 
Constraint ([5]) applied to (jTTJ) implies 



a 2 



^Jv 2 + Aa 2 n(n - 1) - v 
a Q = -— =^> 7(n) = . (12) 



The latter is a convex function with asymptotic behaviour 7(71) "^J 00 ±on and 
symmetric with respect to the quadratic minimum 

n~l/2 yjv 2 — a 2 — v a 2 ( l x 2 



j(n) ~ 1 — , n — - 

V ; 2 V^ 2 - a 2 V 2 

The mean value o.q turns out to be negative but its modulus cannot clearly 
exceed a, which means that only situations with a/v (the analogue of the Kubo 
number) smaller than unity are physically relevant in this picture. Another 
interesting consequence is the relation v\ — v-i = — a, which reduces from three 
to two the number of free parameters in the original definition ([3]) of the strain. 
This means that the lower level of the noise persists longer: in the Lagrangian 
frame, indeed, particles spend longer time in contracting regions than in expand- 
ing ones, thus enhancing the weight of the former in the statistics. Accordingly, 
the Lyapunov exponent, computed as the derivative of (|12[) taken in the origin, 
turns out to be negative [TUl [11] : 

A = 7 '(0) = -a 2 /^. (13) 

The joint probability density function P(R, t) can also be investigated. It 
satisfies the second-order partial differential equation [7 

d 2 P + vd t P - a 2 d R [Rd R {RP)} + a 2 d R (RP) = . (14) 
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Equation (fT¥|) has no stationary solutions, as both P = const, and oc R^ 1 arc 
non-normalizable. It can be solved by means of a saddle-point evaluation after 
variable separation, which amounts to look for the large-deviation asymptotic 
form 

P(ii,i) t_> ~ 00 e- tff W , (15) 

where X = t~ x \n[R/R(0)] can only belong (because of (JTJ) to the interval 
[—a, a]. The Cramer function H(X) is simply the Legendre transform of ~i{n): 



H { X)=»- + X-^^- X2 \ (16) 
V ; 2 2 2a V 7 



It is a convex function vanishing quadratically at the minimum, represented by 
the Lyapunov exponent, H(X) (X — X) 2 v 3 /4a 2 (v 2 — a 2 ), and it approaches 
vertically the boundaries of its compact domain: 



Moreover, it satisfies the Evans-Searles-Gallavotti-Cohen fluctuation relation 
H(X) — H(—X) = X p"2"l[T3] despite the time irreversibility of the strain statis- 
tics [II]. 

Plugging ([16]) into (|15|). we obtain the final resulld 

- [l-V(l-a 2 /^ 2 ){l-ln 2 [-R/i?(0)]/a2 t 2}l vt/2 

P(R,t) t ^°°— 



\R\ 



3 2D case 



In a 2D environment it is possible to consistently impose the incompressibility 
constraint locally, for every single realization of the noise. For the sake of 
simplicity, we shall only consider this situation, which implies the strain matrix 
to be traceless. 

Before studying the general isotropic case, we would like to show a simpler 
instance of flow, namely the hyperbolic one, where a simplification is possible. 

2 It is worth mentioning what would happen if one did not impose the constraint 10 and 
assumed a strain cr(t) = £(t) with the property described in (iii). In this case, 7(n) = 
(V v 2 + 4n 2 a 2 — f)/2 and the Lyapunov exponent would vanish. The Cramer function would 
be H(X) = v/2 - iVa 2 - X 2 /2a, so that 

P(R,t) ' |_R[- 1 exp { - [l - y/l-ln 2 [R/R(0)]/a 2 t 2 Lt/2} 

would be the solution to the equation for the joint probability distribution, similar to H14H 
but without the last term. 
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3.1 Hyperbolic flow 

Let us consider the incompressible irrotational flow described by the diagonal 
strain matrix 

^ -£(t) 

The streamlines are represented by fixed equilateral hyperbolae, which are how- 
ever travelled by fluid particles in randomly-reversing direction. For the sake of 
simplicity the velocity gradients are assumed to have equal reversal probabilities 
and constant absolute values, which justifies the use of the telegraph signal £(£) 
with the properties described in (iii). 

From the analytical point of view, the quantities of our interest which can be 
investigated are of the type 

r„ )fc (i) = R<t(t)R%- k (t) , (17) 

with integer n > k > 0. By means of straightforward algebra, one gets the 
evolution 

I* <+\\ R -\y/v 2 +i{2k-n)*a*+ v ]t/2 , „ J y / v 2 +4(2k-n)2a?-v\t/2 

which implies 

(R n (t)) t_ ~°° e (v / " 2 +4n 2 a 2 -^)t/2 

for positive even n. Therefore, by extrapolation, 



vV 2 + 4n 2 a 2 — v 
7W = j =^ A = , 

which physically reflects the fact that the velocity changes its sign with equal 
probabilities in both directions, thus preventing the interparticle distance from 
exploding exponentially. 

The equation for the joint probability distribution P(R, t) corresponding to (|14j) 
reads 

d 2 t P + vd t P - a 2 d Ri [R\d R] {R)P)\ = , (18) 

with R) = ( 1. J . It is evident that, at long times, such probability must 

factorize into a product of single-coordinate probabilities, in the sense that the 
problem can be reduced to a pair of independent, one-dimensional problems 
described in footnote [H Exploiting this fact, we obtain the final result 

P(R t) e "( 2 -\ /l - ln2 [«i/- R i( )]/ Q2t2 -\/ 1 - ln2 [ fl 2/-R2(o)]/Q 2 t 2 ) I .t/2 n RlR i 
3.2 General isotropic flow 

We now turn to analyse the general isotropic flow described by the strain matrix 

n M _( Ut) Ut) + y/kw(t) \ 
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The noises £i(t), £n(t) and £m(t) are independent of one another but share the 
properties in (iii) with the same coefficients a and v, which can now range from 
to oo independently. Differently from the ID case (see footnoted]), we will see 
that this is enough to mimic a realistic situation. The matrix a thus turns out 
to be singular and nilpotent: deter = = a 2 . 

By means of simple manipulations, analogue to those described for the ID case, 
one can write down a closed system for the coordinates r n j, defined in \Y1\ at 
every n, in the form 



components of r n> . itself, the 3(n + 1) components of r n> . times one noise, the 
3(n + 1) components of r n> . times two (distinct) noises, and finally the n + 1 
components of r n ,- times all the three noises. No other quantity is required 
to close the system, because as soon as a noise appears as square it becomes 
a deterministic quantity and can be taken out of the statistical average. The 
matrix A^"' is thus of order 8(n + 1) x 8(n + 1) and turns out to be dependent 
only on a and v. 

In order to take also rotation into account properly, the presence of three noises 
is necessary and sufficient; one can easily understand that such number would 
sharply increase in higher dimensions, thus strongly enhancing the number of 
components in Rp 1 ^ : this is why we confine ourselves to the 2D case. 
To study the long-time evolution of (r n ,k), one should then look for the eigen- 
values of A( n ) and identify the one with largest positive real part. This is not 
an easy task from the numerical point of view for growing n, and it does not 
look feasible analytically for a generic n, even if the matrix itself can be easily 
written down as a function of n (see appendix [B] for details). In particular, we 
performed the full calculation and we found the whole spectrum for n — 1 and 



The evolution of the linear components and (R2) is very easy to describe: 
the largest eigenvalue is 0, which means that they do not show any exponential 
growth at large times; however, as they can change sign, this does not yield any 
information on (R), i.e. on 7(1). 

The situation is much more meaningful for quadratic components {R\, R1R2, 
i? 2 .), because in this case the largest eigenvalue is given by 




n 



= 2. 



M2 = —v + 



v/216aV + 3V5184aV 2 ~ 3^ 



+ 



\/2lQa 2 v + 3V5184aV 2 ~ 3V 
3 



,1. 



(20) 



which is always positive and behaves asymptotically as 








(21) 
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The eigenvalue [12 expresses the exponential growth rate of (R 2 ) — (Rf) + (i? 2 .) 
and thus coincides with the point 7(2) on the curve 7(71). 

3.2.1 Exact reduced dynamics 

As already pointed out, this general process is computationally very demanding 
and therefore has not been carried out for n > 3. The following step thus con- 
sists in trying to find a reduced dynamics, i.e. a subset of components of TZS 71 ' 
such that the corresponding rows in are nonzero only in columns whose 
index has been taken into account in the reduced subset. 

The first hint into this direction is provided by a more careful look at the matri- 
ces k^ n \ or equivalently at the corresponding system of equations. For instance, 
taking A^ 1 ) into account, one can see (from appendix [B]) that the rows number 1, 
3, 6, 8, 10, 12, 13 and 15 form a closed system, and the same is for the remaining 
8 rows: we have thus reduced from a single 16 x 16 problem to a couple of 8 x 8 
independent problems. The same half-splitting can be easily verified directly 
for A^ 2 ) and is expected to take place also for higher n, thus contributing to a 
strong simplification. 

However, a much heavier simplification can be obtained, even doing without any 
approximation and thus considering exact relations. Here, the key idea is that 
we are not interested in studying the time behaviour in itself of every single com- 
ponent of with l > n + 1 (i.e. of the quantities constructed by multiplying 
the coordinates r n ,fc by some noises), but we simply need some combinations of 
them in order to describe correctly the dynamics of the first n + 1 components; 
moreover, not necessarily all the latter quantities are of our interest, but only a 
combination or a subset of them. 

To be more specific, let us consider again A*- 1 ). Looking at the coefficients which 
appear in the first row, which describes the evolution of Ri , it is natural to con- 
sider the following linear combination: third row, plus sixth row, plus eighth row 
times v2- In other words, we are noticing that the left product of A^ with the 
row vector xo = (1, 0, . . . , 0) gives the vector x\ = (0, 0, 1, 0, 0, 1, 0, \/2, 0, . . . , 0) 
times a. We thus have to perform x\ ■ , and it can easily be seen that the result 
is —vx\. Basically, we are down to a 2 x 2 matrix (note that 2 = l + l = n+l) 
with zeros on one line and (a, —v) on the other: its eigenvalues are trivially and 
—v\ being the larger one, in accordance with the result from the complete form. 
It is clear that the same result also holds considering the evolution of Ri (second 
line of AW), i.e. x = (0, 1, 0, ... , 0) and x t = (0, 0, 0, -1, 1, 0, -y/2, 0, . . . , 0). 
Let us turn to A^ 2 ). Our aim is to find the behaviour of R 2 , therefore let us begin 
with adding up the first and third lines, corresponding to R\ and R\ respec- 
tively. In formulas, if yo = (1, 0, 1, 0, . . . , 0), we have yo • A^ 2 ' = 2ay\, with y\ = 
(0, 0, 0, 1, 0, -1, 0, 2, 0, . . . , 0). Then we compute yi ■ A (2 *> = Aay - vyi + 2V2ay 2 , 
where j/2 = (0, . . . , 0, 2, 0,-1, 0, 1, 0, 0, 0). Finally, we are able to close the sys- 
tem with 2/2 • A^ 2 ' = — 2\[2ay\ — Ivy^- The resulting 3x3 matrix (notice that 
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3 = 2 + l= n + l)is thus 



2a \ 

4a -v 2-y/2a I 
-2V2a -2i/ / 

and its eigenvalues satisfy 

H% + Zvnl + 2v 2 H2 - I6a 2 v = , (22) 

which gives the value of H2 reported in (|2H)) as the solution with largest positive 
real part. It is also worth noticing that the highest nonzero component in the 
vectors i/OjZ/i>2/2 is the last but three, which means that the final 3 components 
of VS 2 ^ (which correspond to the coordinates r^^k times all three noises) are not 
involved at all in the dynamics of (R 2 )- 

Having tested this simplification to (n + 1) x (n+1) on the already-solved cases 
n = 1 and n = 2, we now turn to higher n, but we confine ourselves to even 
values. Indeed, 



R n = (Rf + R\) nl2 = 




by definition for positive even n, thus the corresponding points 7(n) can be 
extracted from the knowledge of n/2 + 1 of the first n+1 components in Tv- n '. 
On the contrary, this is not the case for the odd n's, which behave similarly to 
non-natural n's because they involve nonlinear operations (in this case a square 
root) that do not commute with the statistical average (•). 
In particular, we analysed the values n = 4, n = 6 and n = 8. Adopting a tech- 
nique similar to the one described in the previous paragraphs, one gets matrices 
of order 5, 7 and 9 respectively (i.e. n + 1), whose characteristic polynomials 
are: 

l4 + 8^4 + 23t^| + (28^ 3 - 96a 2 v)ni + {l2v A - 336(rV 
+ (-288a 2 ^ 3 - 1152aV) , 

(4 + Uvf4 + A9v 2 ^l + (113z/ 3 - 96a 2 v)(ij + (142*/ 1 - 768a 2 i/ 2 )/4 
+ (92v 5 - 2208a 2 ^ 3 - 13440a'V)^ 2 ; + (24i/ 6 - 2688a 2 ;/ 1 - 44160a V)/i 6 
+ (-1152a 2 ^ 5 - 34560a 4 !/ 3 - 230400a e V) 

and 

f4 + 12v(4 + 60^ 2 u| + (128a 2 ^ + 162i/ 3 )/^ + (255^ 4 + 960aV)/xf 
+(234u 5 + 2816a 2 ^ 3 - 44160a 4 ^)^ + (116v 6 + 4032a 2 */ 1 - 256384aV)u 3 
+{2Av 7 + 2816aV - 567552a 4 ^ 3 - 3548160aV)/i 2 ! + (768a 2 i^ 6 - 570368a 4 i/ 
-10407936a f V 2 )^ 8 + (-215040aV - 7454720a f V 3 - 92897280a 8 ^) . 
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The respective roots with largest positive real part behave asymptotically as: 



Hi 



A*6 



^8 



24 576— for - 

V V V 

for - 
v 



48^-2016^ for- 
a 



oo 



2v / 1800a 6 i/ for - 
v 

80— - 5280^ for- 

V V V 

a 
v 



2v / 181440a 8 !/ for 



These relations provide us with three points of the curve 7(n) in the asymptotic 
conditions, to which three more can be added: besides the already-mentioned 
7(2), one should indeed remember that both 7(0) and 7(— 2) (in 2D) must 
vanish [T5J, Unj [T] . A simple extrapolation then suggests the following asymptotic 
behaviours: 

G(n)-+0(n)4 for-^0 
7(n) ~ { v v v (24) 

g( n )a n K n+r >v 1 K n+1 ') for - -> 00 . 

v 

The coefficient Gin) can be proven rigorously to be n(n + 2) (see the next sub- 
section), while such a proof does not look feasible for Q(n). However, starting 
from the quadratic form of the former, one can guess a fourth-order polynomial 
for the latter, G{ri) = 3n(n + 2)(n 2 + 2n + 8)/4, and the correctness of this 
expression seems to be confirmed by the fact that it satisfies the six conditions 
in n = —2,0,2,4,6,8 despite possessing only five degrees of freedom. As the 
following coefficients in the power expansion at small ajv are expected to show 
higher powers of n, such an expansion is not uniform in the sense that it must 
fail for some (large enough) value of n. On the contrary, the coefficient g(n) 
must originate from some combinatorics, but has not been identified yet, while 
the presence of non-integer powers of a and v suggests that the large-a/j/ devel- 
opment is singular. It is worth noticing the dependence of the growth rate of 
the n-th moment as a function of the noise correlation time (inverse of v) : we 
have a linear growth when the latter is small and a power-law decay when it is 
large. 

From (j24|) . exploiting the knowledge of G(n) and the vanishing of (?(0), we get 
the asymptotic behaviours of the (positive, as expected) Lyapunov exponent: 

2 

a a 
2— for > 

v v (25) 

a 

g'(0)v for ► 00 . 

v 
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3.2.2 Quasi delta-correlated case 

The results obtained in the previous subsection are exact, but they have been 
obtained for specific values of n, and the generalization of the coefficients (|24j) 
as functions of n was empirical. There is however at least one instance in which 
rigorous analytical results can be obtained: it is provided by the quasi delta- 
correlated case, corresponding to small a/ v. 

It is worth reminding that in the exactly delta-correlated case [TJ [5] the prob- 
ability distribution of R(t) is lognormal, the Cramer function is quadratic and 
the Lyapunov exponent is given by 2Di, where D\ is the multiplicative constant 
in the well-known expression for the spatial part of the two-point structure func- 
tion of the velocity field. A comparison with (|25|) suggests the correspondence 
a 2 jv «-> D%, which means that the correct limit is a — > oo <— v such that 
ajv — » but a? jv keeps finite. 

The key point in this limit is that the leading behaviour can be extracted rig- 
orously, for any natural n, from a restricted dynamics, namely keeping into 
account only the first 4(n + 1) components in TZS n > (corresponding to the aver- 
age of the coordinates r n ,k, an d of the latter times one noise). To show this, for 
any quantity Y let us define its small-Kubo expansion as 

oo , 
(=0 

and let us plug it in the evolution equation (fT9|) . 
At the lowest order (1 = 0), we obtain 

for 1 < l < n + 1 
for (n+l) + 1 < b < 4(71 + 1) 
for 4(n+ 1) + 1 < l < 7(n+ 1) ^ 
for 7(n + l) + 1 < i < 8(n + l) , 

so that (Tl^i 1 ) is a constant in the group of its first n + 1 components (corre- 
sponding to the average of the coordinates r n> fc) and decreases exponentially in 
the remaining three groups. 

The equations for I = 1 have the same operatorial structure as in (p?T|) but, 
instead of being homogeneous, are rather forced by the aforementioned I = 
quantities in a "cross fashion", i.e. not by quantities of the same group but 
of "nearest-neighbour" groups. Therefore, the only nondecreasing source term 
(C^roil) with 1 < t < n + 1) appears exclusively in the evolution equations for 
with (n+l) + l<K< 4(n +1), which implies that now (TZ^) evolves 
towards a finite constant in its first 4(n + 1) components (corresponding to the 
average of the coordinates r n ,k, and of the latter times one noise) and decreases 
exponentially in the remaining ones. 

Moving to / = 2, the first 7(n + 1) components are now forced by a non- 
vanishing ("constant") source term. This means that ) goes to a constant 
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for (n+ 1) + 1 < l < 7(n + 1), while the first n+ 1 components (corresponding 
to the average of the coordinates r^) get a linear behaviour in time. For / = 3 
also the components with (n+1) + 1 < l < 4(n + 1) behave linearly, which in 
turn implies a quadratic behaviour for the first n+1 components at I = 4, and 
so on. 

The situation is resumed in the following table, in which for every I we show the 
degree T of the leading behaviour in time of each group of components 
E.g., denotes the evolution towards a finite constant, 1 a linear behaviour and 
x the only presence of a decreasing exponential. 



1 = 





1 


2 


3 


4 


5 


6 


1 < i < n + 1 










1 


1 


2 


2 


3 ... 


i-e. (r„,fc) 


















(n+ 1) + 1 < t < 4(rc- 


f 1) 


X 








1 


1 


2 


2 ... 


i.e. (£.r njfe ) 


















4(n+l) + l < i < 7(n 
i.e. (£,£, r n , k ) 


+ 1) 


X 


X 








1 


1 


2 ... 


7(n+l) + l< i < 8(n 
i.e. (£i£ii£iiirn,k) 


+ 1) 


X 


X 


X 








1 


1 ... 



Analysing the first line of this table, we argue that the expected exponential 
behaviour for {r n ,k) corresponds to the summation of the series in (f2"o| in the 
fashion 



00 2L 

^W) = E(3 (vt) L f(n,k,L) [l + 0Q] '^V"* (for some fc) 



(28) 

One would be tempted to conclude that, for each n, the largest-positive-real- 
part eigenvalue fj, n could be easily extracted by looking at the coefficient in (|28[) 
with L = 1, corresponding to the column Z = 2 in the table, because e ce ~ 1 + ce 
for small e. However, this would be the case only if the eigenvalues correspond- 
ing to a fixed n were definitely separate for ajv — > 0. On the contrary, from 
the study of the full matrix performed e.g. for n = 2 (see also appendix [Bj, we 
know that several eigenvalues vanish in this limit, and we have to identify the 
one which keeps the largest positive real part in this process: we are therefore 
facing a degenerate perturbation theory, and the quest for the asymptotic ex- 
pression of fi n (or, in a way, for the rigorous proof that it has the expression 
shown in the previous subsection) is somehow longer. 

Let us then proceed as follows. From the inspection of the system (|27l) and of 
the following orders, it is apparent that the structure of the equations is such 
that the degree T of the leading behaviour in time (i.e. the quantity reported in 
the table) can grow of one unity only when a quantity of the second group, of 
order I and time degree T, acts as a source term for a quantity of the first group 
of order l + l, which therefore gets a degree T + 1. This means that, in order to 
find the leading behaviour of \x n at small a/v, one can completely neglect the 
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lower two rows in the table and only consider the upper two. Indeed, starting 
from the only non-vanishing group (T = 0) in the column / = and remem- 
bering that the forcing mechanism acts in a diagonal "cross fashion" (i.e. every 
term can only force its upper-right and lower- right neighbours), the components 
of the third and of the fourth group can give rise to linear terms in the first 
group only for I = 4 and / = 6 respectively, which are clearly subdominant with 
respect both to the quadratic and cubic behaviours of these same components 
(T = 2 and T = 3 in the table) and to the linear behaviour found for I = 2 due 
to the forcing by the second group. Back to the initial dynamics d t (TZ^) = ■ ■ ., 
this amounts to say that in the evolution equation for {r n> k) one should keep the 
complete form into account, while for {(,»r n ^) one can neglect all source terms 
involving two different noises: the system thus closes at this stage, reducing the 
order of the corresponding matrix from the initial 8(n + 1) to A{n + 1). In other 
words, only the upper-left quarter of the matrix A*™) is relevant. 
The problem can now easily be recast as a system of n + 1 second-order differ- 
ential equations for r n ^. Indeed, we have 

d 2 r nM = kR k - 1 R^ k d 2 R 1 + (n - k)R k R%- k - 1 d 2 R 2 

+k(k - ^R^R^idtRi) 2 + (n-k){n-k- l)R k R%- k - 2 (d t R 2 ) 2 
+2k(n- k)R k - 1 R 2 l - k - 1 d t Rid t R 2 , 

and, keeping into account the aforementioned simplification in the equations for 

d 2 (r n , k ) + vd t (r n>k ) = a 2 [(n 2 + 6k 2 - 6nk - n)(r n , fc ) (29) 
+3k(k - l)(r n , fe _ 2 ) + 3(n- fe)(n - k - l)(r„^ +2 )] 

= a2 Mfcfc' ( r n,k') ■ 

The matrix M^ n \ defined by the right-hand side of ([29]) . is a priori of order 
(n +1) x (n + 1). However, it is apparent from (f29|) that, for a fixed n, the 
dynamics of the components of {r n ,k) with even k are independent of those with 
odd k. Therefore, as our main goal is to reconstruct objects like (|23|) . we can just 
focus on subset IVH n ) of M^") consisting of the lines with even k, and thus reduce 
to a (n/2 + 1) x (n/2 + 1) problem, simply neglecting the lines corresponding 
to odd k. In other words, we rewrite (f2"9"| as 

d 2 t {r nak )+vd t {r n , 2k )=a 2 £ M«, (r„, fc /) = a 2 u[ n k ] (r n>w ) . (30) 

A:' =0,2 n 

Moreover, interpreting (|30|) as 

{d t (?V2fc) = -^(r„,2fe) + vY n £k 
d t Y nak = ^Mj$(r n , 2 fc<> , 

where Y n> 2k = (fn,2k) + dt(r ni2k ) /v, it is easy to show that the largest-positive- 
real-part eigenvalue /i„ (divided by a 2 jv) that we are looking for must belong to 



13 



the spectrum of the subset matrix MM, because the remaining eigenvalues are 
all — v. Such subset matrix is tridiagonal centrosymmetric but not symmetric, 
which means that left and right eigenvectors differ. It is a simple task to prove 
that the value n(n + 2) (guessed for // 2 /(a 2 /^) asymptotically from the previous 
subsection) is actually an eigenvalue, but it is not possible to use a variational 
method to show that it has the largest positive real part. However, this is not 
needed, because the row vector built with the n/2+1 components of the binomial 
coefficient ( n [ 2 ), Vk = (1, n/2, . . . , n/2, 1) for k = 0, 1, . . . , n/2, happens to be 
the left eigenvector of M^™) corresponding to the eigenvalue n(n + 2). Therefore, 
using (|23|) and ([30]) , one has: 

E ^(d 2 t(r n ,k)+vd t (r ntk )) 



d 2 t (R n ) +vd t {R n 



fc=0.2,. 



E 



n l j (d 2 (r n>2 k) + vdt{r n ,2k)) 



fc=0,l,...,n/2 
fe=0,l,...,n/2 fe'=0,l,...,n/2 



This implies 



-« 2 E E ^ \ (r n , 2k ,) 

fc'=0,l,...,n/2 \k=0,l,...,n/2 

= a 2 22 n(n + 2)v k > (r ni 2k>) 
fc'=0,l,...,n/2 

f nj2^ 
k' 



n(n + 2)a 2 ^ 

fe'=0,l,...,n/2 



\r n ,2k' 



= n(n + 2)a 2 V] 



fc'=0,2,...,n 



n/2 
k'/2 



[r n ,k> 



= n(n + 2)a 2 (R n ) 



Mi 



a/i/->0 ^/v 2 + 4n(n + 2)a 2 - v 



^2 r 



= n(n + 2)- 



1 + 



as expected. 

Notice that the next-leading correction cannot be captured at this stage, because 
of the approximation introduced previously. To take it into account correctly, 
one should reformulate the analysis performed in this section, including also 
quantities like (£»£orVi,fc) but however excluding (£,i£,n£,nir n , k ), i.e. considering 
a reduced dynamics in the first 7(n + 1) components of lZS n ' . 
It is worth mentioning that an analysis like the previous one, but applied to odd 
n's, leads to 

2 

a/V->0 . , , ,a 

u„ ~ (n - l)(n + 3) — , 
v 
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in accordance e.g. with the exact result \i\ = found with the complete dynam- 
ics. However, as already pointed out, such values are not related to the curve 
7(71). 

Lastly, one would be interested in reformulating the previous study in the 
quasi-deterministic limit, a/v — > 00 (see appendix Unfortunately, this is 
not possible, because the limit turns out to be singular: an exponential growth 
is expected in general for finite a/v, but a power law is found when this ratio 
is infinite. In other words, an expansion like (|26p with a/v inverted does not 
work. In this reflects into the presence of non-integer powers and into the 
difference between the limits a — > 00 and v — ► 0. 

4 Conclusions 

By assuming the telegraph-noise model for the velocity gradient (or strain ma- 
trix), we were able to carry out analytical computations and to obtain sev- 
eral results on the separation between two fluid particles. Focusing on smooth 
flows (Batchelor regime), we firstly analysed the one-dimensional compressible 
case, finding explicit expressions for the long-time evolution of the interparticle- 
distance moments fj 12[) . the Lyapunov exponent (113j) and the Cramer function 
(|16p . Then wc concentrated on the two-dimensional incompressible situation 
and provided an elementary extension of the ID case, represented by the hy- 
perbolic flow. Moving to the general isotropic thorough analysis on the 
complete dynamics was made for the evolution of linear (n = 1) and quadratic 
(n = 2) components. Due to high computational cost, at higher n we focused on 
a restricted, though exact, dynamics: however, only specific values of n could be 
studied in this way, leading to the extrapolations (pM|) and (j2"5| . Such guess was 
rigorously proved in the quasi-delta-correlated limit, for which approximated 
equations were introduced. 

We believe that the present paper represents an interesting example of the 
use of a coloured noise for which the well-known closure problem can be solved 
analytically [17] , 
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A Appendix: deterministic case in 2D 

The deterministic case corresponds to v = 0, when the noise takes a constant 
value ±o. In 2D, the presence or absence of rotation plays a crucial role. Indeed, 
in the diagonal case (corresponding to the hyperbolic flow shown in section l3"7Tj) , 

one finds the exponential evolution R n e ant , which implies 7(71) = an and 
A = a. 

On the contrary, in the general isotropic case (with rotation), one gets d 2 Ri — 
= d^ i?2 (for any of the 2 3 = 8 possible combinations of the noise signs), thus 
both components are linear in time and R n t n . Consequently, a power- 
law temporal dependence is found for the separation and A = 0. It is worth 
noticing that this result is due to the fact that, with our choice, a 2 = and is 
characteristic of the 2D isotropic situation. 
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B Appendix: calculation details for 2D 

In this appendix we provide some details of the calculation for the 2D general 

isotropic case. 

The matrix k$ reads: 

[(n + 1) - 2i + l]a^ K , t +(„+i) + [(n + 1) - t]a<5 K , t+2 („+i)+i 
+V2[(n + 1) - t]a5 K)( , + 3( n+ i) + i + (a - l)aS K>l+2 (n+i)-i 
-V2(i - l)aS K ^ +3(n+1) _ 1 

for 1 < i < n + 1 , 

[3(n + 1) - 2a + l]a£ Kjt _(„ +1) + [2(n + 1) - t]a<5 Kit+3( „ + i )+ i 
+V2[2(n + 1) - i]a5 K . L+A(n+1)+1 + [-(n + 1) + a - l]a£ K)t+3 ( n+ i)_i 
-V2[-(n + 1) + a - l]a5 Kil+4 („ + i)_i - z/£ K , t 

for (n+ 1) + 1 < a < 2(n+ 1) , 

[5(n + 1) - 2a + l]o5 K)t+ 2(„+i) + [3(n + 1) - t]a5 K)t _ 2 ( n+ i)+i 
+V2[3(n + 1) - A]a<5 K , t+4(rl+1)+1 + [-2(n + 1) + a - l]a(5 Kjt _ 2 („+i)-i 
-V2[-2(n + 1) + a - l]a<S re)( ,+4(n+i)-i - 

for 2(n+ 1) + 1 < a < 3(n+ 1) , 

[7(n + 1) - 2a + l]a5 K)t+ 2(„+i) + [4(n + 1) - t]a£ Kit+3( „ + i) + i 
+V2[4(n + 1) - t]a5 K)t _ 3 ( n+ i) + i + [-3(n + 1) + i - l]a < 5 K ^ +3(rl+1) _ 1 
-V2[-3(n + 1) + a - l]a(5 Kit _ 3( „ + i)_i - v5 K<L 

for 3(n+ 1) + 1 < a < 4(n+ 1) , 

< 

[9(n + 1) - 2a + l]a^ it _ 2 („+i) + [5(n + 1) - A]a5 Kit _ 3(n+1)+1 
+V2[5(n + 1) - t]a5 K)t+3(n+ i )+ i + [-4(n + 1) + a - l]^,,.^.^.! 
-V2[-4(n + 1) + a - l]a(5 Kit+3 („ + i)_i - 2i/<5 K , t 

for 4(n+ 1) + 1 < a < 5(n+ 1) , 

[ll(n + 1) - 2a + l]a5 K)t _2(„+i) + [6(n + 1) - t]o5 K)( , + 2( n +i)+i 
+V2[6(n + 1) - t]a5 K)t _ 4 ( n+ i) + i + [-5(n + 1) + a - l]a£ KiH _ 2 („+i)_i 
-V2[-5(n + 1) + a - l]a<5 Kit _ 4 („ + i)_i - 2i/<S K)( , 

for 5(n+ 1) + 1 < a < 6(n+ 1) , 

[13(n + 1) - 2a + l]a£ K)t+(n+1 ) + [7(n + 1) - i]a£ K)t _ 3(n+1 ) +1 
+V2[7(n + 1) - i]a£ K)t -4(n+i)+i + [-6(n + 1) + i - l]a£ K , t - 3 („+i)-i 
— v/2[-6(n + 1) + a - l]a5 Kit _ 4 („ + i)_i - 2^<5 K , t 

for 6(n+ 1) + 1 < a < 7(n+ 1) , 

[15(n + 1) - 2a + l]a^ jt _ ( „ + i) + [8(n + 1) - i]a5 K ^_ 2(n+1)+1 
+V2[8(n + 1) - t]a5 K)t _ 3 ( n+ i) + i + [-7(n + 1) + i - l]a5 K)t _ 2 ( n +i)-i 
-\/2[-7(n + 1) + a - l]a5 re , t _ 3 („ + i)_i - 2^<5 K , t 

for 7(n+ 1) + 1 < a < 8(n+ 1) . 

We focus at first on the study of linear coordinates (n = 1): we introduce 
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the 16-component vector 



/ Ri \ 

R2 
£i(t)Ri/a 
^(t)R 2 /a 
£n(t)Ri/a 
£n{t)R 2 /a 
Zin(t)Ri/a 
= t,in(t)R 2 /a 

Ut)Ut)Ri/a 2 
Ut)tn(t)R2/a 2 

Ut)^ii(t)Ri/a 2 
ti(t)tiiii(t)R2/a 2 
Zii{t)Ziii{t)Ri/a 2 
Zn(t)!iiii(t)R 2 /a 2 
Ut)tnmiii(t)Ri/a 3 
V Ut)Znmm(t)R2/a 3 J 

The matrix A^ 1 ^ reads 
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a 





V2a 
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-V2a 
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V2a 
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— v 
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— V 
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a 














— V 











—a 
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-V2q 
















V2a 














— V 











a 








a 










V2a 




















— V 











—a 


a 
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a 











-2v 




















V2a 










a 








—a 











-1v 
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V2a 








a 











-2v 














a 
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V2a 
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a 






















V2a 





a 














-1v 





a 
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-V2a 
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-2i/ 





—a 































V2a 
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-Zv 





V 
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-V2a 





a 








—a 





-3i7 



and its eigenvalues are given by (twice), — v (six times), —2v (six times) and 
—iv (twice). Therefore, the largest eigenvalue is fix — 0, with eigenvectors 
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via 
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v a 
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-1 















1 
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V2 






-V2 
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Let us now move to quadratic quantities (n — 2): we define the 24-component 
vector 

/ Rl \ 

R\R 2 

Ut)RiR 2 /a 

Zn{t)R 2 /a 
Zn{t)RiRa/a 
t,n{t)Rl/a 
ini(t)R\/a 
$m(t)R 1 R 2 /a 
7? (2) _ £,in{t)Rl/a 

Ut)Zn(t)Ri/a 2 
Ut)Zii{t)RiR2la 2 

Umm(t)Ri/a 2 
Ut)t;iii(t)RiR2/a 2 

Ut)Zin(t)R 2 2/a 2 

Zii{t)iin{t)R 2 /a 2 
Cii(t)On(t)RiR 2 /a 2 

in{t)ini{t)Rl/a 2 

umnmin(t)Ri/a 3 

Ut)Zii(t)Ziii(t)RiR2/a 3 

V Ut)£,ii{t)iin{t)Rll^ J 

The associated matrix is 
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and its eigenvalues are: —2v (three times), — v (three times), the three roots of 
//I + 3i^yu| + 2v 2 [i2 — lQa 2 v = (corresponding to the reduced dynamics (|22)) ). 
the three roots (each taken twice) of + 3^/x 2 , + 2v 2 [ii + Sa 2 v — 0, the three 
roots (each taken twice) of /z?, + + llv 2 ^ + (6^ 3 — 8a 2 ^) = and the three 
roots of ^2 + 6vfi 2 + Wv 2 [i2 + (6^ 3 + 16a 2 v) — 0. In accordance with (|20|) . the 
largest eigenvalue is 

v 2 ^2lQa 2 v + 3V5184a 4 ^ 2 

H2 = -v+ ~, H „ 

V 21&a 2 v + 3V5184a 4 ^ 2 - 3z/ 6 3 

and its eigenvector reads: 

/ n/2 + 2\/§i///i \ 


V~2 + 2\/2v/n 
v/V2a + /i/2\/2a 


-u/V2a- n/2\/2a 


^/\/2a + /V 2 V2a 
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